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This paper proposes rigorous geometrical treatment of bioelectrodynamics, under- 
pinning two fast-growing biomedical research fields: hioelectromagnetism, which deals 
with the ability of life to produce its own electromagnetism, and bioelectromagnetics, 
which deals with the effect on life from external electromagnetism. 
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1 Introduction 

Bioelectromagnetism, sometimes equated with bioelectricity, bioelectromagnetism refers to 
the electrical, magnetic or electromagnetic fields produced by living cells, tissues or or- 
ganisms. Examples include the cell potential of cell membranes and the electric currents 
that flow in nerves and muscles, as a result of action potentials. It is studied primar- 
ily through the techniques of electrophysiology (see, e.g. p]), which is the study of the 
electrical properties of biological cells and tissues. It involves measurements of voltage 
change or electrical current flow on a wide variety of scales from single ion channel pro- 
teins to whole tissues like the heart. In neuroscience, it includes measurements of the 
electrical activity of neurons, and particularly neural action potentials, or neural spikes, 
pioneered by five celebrated Hodgkin-Huxley papers in 1952 [U [3l IU (Sj [6] (formalized in 
[6] by a 1963-Nobel-Prize winning mathematical model), which are pulse-like waves of 
voltage that travel along several types of cell membranes. The best-understood example 
is generated on the membrane of the axon of a neuron, but also appears in other types 
of excitable cells, such as cardiac muscle cells, and even plant cells|l| The principal ions 



^The resting voltage across the axonal membrane is typically -70 millivolts (mV), with the inside being 
more negative than the outside. As an action potential passes through a point, this voltage rises to roughly 
+40 mV in one millisecond, then returns to -70 mV. The action potential moves rapidly down the axon, 
with a conduction velocity as high as 100 meters/second (224 miles per hour). Because they are able to 
transmit information so fast, the flow of action potentials is a very efficient form of data transmission, 
considering that each neuron the signal passes through can be up to a meter in length. An action potential 
is provoked on a patch of membrane when the membrane is depolarized sufficiently, i.e., when the voltage 
of the cell's interior relative to the cell's exterior is raised above a threshold. Such a depolarization opens 
voltage-sensitive channels, which allows current to flow into the axon, further depolarizing the membrane. 
This will cause the membrane to 'flre', initiating a positive feedback loop that suddenly and rapidly causes 
the voltage inside the axon to become more positive. After this rapid rise, the membrane voltage is restored 



involved in an action potential are sodium and potassium cations; sodium ions enter the 
cell, and potassium ions leave, restoring equilibrium. Relatively few ions need to cross the 
membrane for the membrane voltage to change drastically. The ions exchanged during 
an action potential, therefore, make a negligible change in the interior and exterior ionic 
concentrations. The few ions that do cross are pumped out again by the continual action 
of the sodiumpotassium pump, which, with other ion transporters, maintains the normal 
ratio of ion concentrations across the membraneo Calcium cations and chloride anions 
are involved in a few types of action potentials, such as the cardiac action potential |9j. 



to its resting value by a combination of effects: the channels responsible for the initial inward current are 
inactivated, while the raised voltage opens other voltage-sensitive channels that allow a compensating 
outward current. Because of the positive feedback, an action potential is all-or-none; there are no partial 
action potentials. In neurons, a typical action potential lasts for just a few thousandths of a second at 
any given point along their length. The passage of an action potential can leave the ion channels in a 
non-equilibrium state, making them more difficult to open, and thus inhibiting another action potential at 
the same spot: such an axon is said to be refractory. 

^The action potential 'travels' along the axon without fading out because the signal is regenerated at 
each patch of membrane. This happens because an action potential at one patch raises the voltage at 
nearby patches, depolarizing them and provoking a new action potential there. In unmyelinated neurons, 
the patches are adjacent, but in myelinated neurons, the action potential 'hops' between distant patches, 
making the process both faster and more efficient. The axons of neurons generally branch, and an action 
potential often travels along both forks from a branch point. The action potential stops at the end of 
these branches, but usually causes the secretion of neurotransmitters at the synapses that are found there. 
These neurotransmitters bind to receptors on adjacent cells. These receptors are themselves ion channels, 
althoughin contrast to the axonal channelsthey are generally opened by the presence of a neurotransmitter, 
rather than by changes in voltage. The opening of these receptor channels can help to depolarize the 
membrane of the new cell (an excitatory channel) or work against its depolarization (an inhibitory channel). 
If these depolarizations are sufficiently strong, they can provoke another action potential in the new cell. 
The flow of currents within an axon can be described quantitatively by cable theory [7] (as demonstrated 
by Hodgkin [8]). In simple cable theory, the neuron is treated as an electrically passive, perfectly cylindrical 
transmission cable, which can be described by a PDE: 

dt dx'^ 

where V{x, t) is the voltage across the membrane at a time t and a position x along the length of the 
neuron, and where A and r are the characteristic length and time scales on which those voltages decay in 
response to a stimulus. 

^The cardiac action potential differs from the neuronal action potential by having an extended plateau, 
in which the membrane is held at a high voltage for a few hundred milliseconds prior to being repolarized 
by the potassium current as usual. This plateau is due to the action of slower calcium channels opening 
and holding the membrane voltage near their equilibrium potential even after the sodium channels have 
inactivated. The cardiac action potential plays an important role in coordinating the contraction of the 
heart. The cardiac cells of the sinoatrial node provide the pacemaker potential that synchronizes the heart. 
The action potentials of those cells propagate to and through the atrioventricular node (AV node), which 
is normally the only conduction pathway between the atria and the ventricles. Action potentials from the 
AV node travel through the bundle of His and thence to the Purkinje fibers. Conversely, anomalies in the 
cardiac action potential, whether due to a congenital mutation or injury, can lead to human pathologies. 



The action potential in a normal skeletal muscle cell is similar to the action potential in 
neurons |271 [22] . Action potentials result from the depolarization of the cell membrane (the 
sarcolemma), which opens voltage-sensitive sodium channels; these becomes inactivated 
and the membrane is repolarized through the outward current of potassium ions. The 
resting potential prior to the action potential is typically -90mV, somewhat more negative 
than typical neurons. The muscle action potential lasts roughly 24 ms, the absolute 
refractory period is roughly 13 ms, and the conduction velocity along the muscle is roughly 
5 m/s. The action potential releases calcium ions that free up the tropomyosin and allow 
the muscle to contract. Muscle action potentials are provoked by the arrival of a pre- 
synaptic neuronal action potential at the neuromuscular junction. There output is muscle 
force generation. 

Bioelectromagnetic activity of the brain is measured by non-invasive electroencephalog- 
raphy (EEG), usually recorded from electrodes placed on the scalpj (see, e.g. [TT1) . Nor- 
mal brain wave patterns are classified by the frequency range as: delta (< 4Hz)j^ theta 
(4 - 7Hz)|l alpha (7 - 12Hz) III beta (12 - 30Hz)lf and gamma (30 >Hz)|!| 

Bioelectromagnetic activity of the heart is measured by non-invasive electrocardiog- 



especially arrhythmias (see, e.g. [lOj). 

''The data measured by the scalp EEG are used for chnical and research purposes. A technique similar to 
the EEG is intracranial EEG (icEEG), also referred to as subdural EEG (sdEEG) and electrocorticography 
(ECoG). These terms refer to the recording of activity from the surface of the brain (rather than the scalp). 
Because of the filtering characteristics of the skull and scalp, icEEG activity has a much higher spatial 
resolution than surface EEG. 

^ Delta is the highest in amplitude and the slowest in frequency. It is seen normally in adults in slow 
wave sleep. It is also seen normally in babies. It may occur locally with subcortical lesions and in general 
distribution with diffuse lesions, metabolic encephalopathy hydrocephalus or deep midline lesions. It is 
usually most prominent frontally in adults (e.g. Frontal Intermittent Rhythmic Delta) and posteriorly in 
children e.g. Occipital Intermittent Rhythmic Delta). 

^ Theta is seen normally in young children. It may be seen in drowsiness or arousal in older children 
and adults, as well as in meditation. 

'^Alpha is seen in the posterior regions of the head on both sides, being higher in amplitude on the 
dominant side. It is brought out by closing the eyes and by relaxation. It was noted to attenuate with 
eye opening or mental exertion. This activity is now referred to as 'posterior basic rhythm', the 'posterior 
dominant rhythm' or the 'posterior alpha rhythm', and is actually slower than 8Hz in young children 
(therefore technically in the theta range). In addition, there are two other normal alpha rhythms that are 
typically discussed: the Mu rhythm (which is an alpha-range activity that is seen over the sensorimotor 
cortex and characteristically attenuates with movement of the contralateral arm) and a temporal 'third 
rhythm'. 

^Beta is usually seen on both sides of the scalp in symmetrical distribution and is most evident frontally. 
Low amplitude beta with multiple and varying frequencies is often associated with active, busy or anxious 
thinking and active concentration. 

^Because of the filtering properties of the skull and scalp, gamma rhythms can only be recorded from 
electrocorticography or possibly with magnetoencephalography. Gamma rhythms are thought to represent 
binding of different populations of neurons together into a network for the purpose of carrying out a certain 
cognitive or motor function. 



raphy (ECG or EKG), usually recorded from electrodes placed on the skin of a thorax 
(see, e.g. ^2J)1^ Human heartbeats originate from the sinoatrial node (SA node) near 
the right atrium of the heart. Modified muscle cells contract, sending a signal to other 
muscle cells in the heart to contract. The signal spreads to the atrioventricular node (AV 
node). Signals carried from the AV node, slightly delayed, through bundle of His fibers and 
Purkinje fibers cause the ventricles to contract simultaneously. ECG measures changes in 
electrical potential across the heart, and can detect the contraction pulses that pass over 
the surface of the heart. There are three slow, negative changes, known as P, R, and T. 
Positive deflections are the Q and S waves. The P wave represents the contraction impulse 
of the atria, the T wave the ventricular contraction!^ 



^"Electrodes on different sides of the heart measure the activity of different parts of the heart muscle. 
An ECG displays the voltage between pairs of these electrodes, and the muscle activity that they measure, 
from different directions, also understood as vectors. This display indicates the overall rhythm of the 
heart, and weaknesses in different parts of the heart muscle, ft is the best way to measure and diagnose 
abnormal rhythms of the heart, particularly abnormal rhythms caused by damage to the conductive tissue 
that carries electrical signals, or abnormal rhythms caused by levels of dissolved salts (electrolytes), such as 
potassium, that are too high or low. In myocardial infarction (MI), the ECG can identify damaged heart 
muscle. But it can only identify damage to muscle in certain areas, so it can't rule out damage in other 
areas. The ECG cannot reliably measure the pumping ability of the heart; for which ultrasound-based 
(echocardiography) or nuclear medicine tests are used. 

^^The baseline voltage of the electrocardiogram is known as the isoelectric line. A typical ECG tracing 
of a normal heartbeat (or, cardiac cycle) consists of a P wave, a QRS complex and a T wave (see, e.g. 
[13)). A small U wave is normally visible in 50 to 75% of ECGs. During normal atrial depolarization, 
the mean electrical vector is directed from the SA node towards the AV node, and spreads from the right 
atrium to the left atrium. This turns into the P wave on the EKG, which is upright in II, III, and aVF 
(since the general electrical activity is going toward the positive electrode in those leads), and inverted in 
aVR (since it is going away from the positive electrode for that lead). The relationship between P waves 
and QRS complexes helps distinguish various cardiac arrhythmias. The shape and duration of the P waves 
may indicate atrial enlargement. The PR interval is measured from the beginning of the P wave to the 
beginning of the QRS complex. It is usually 0.12 to 0.20 s (120 to 200 ms). A prolonged PR interval may 
indicate a first degree heart block. A short PR interval may indicate an accessory pathway that leads to 
early activation of the ventricles, such as seen in Wolff-Parkinson- White syndrome. A variable PR interval 
may indicate other types of heart block. PR segment depression may indicate atrial injury or pericarditis. 
The QRS complex is a structure on the ECG that corresponds to the depolarization of the ventricles. 
Because the ventricles contain more muscle mass than the atria, the QRS complex is larger than the P 
wave. In addition, because the His/Purkinje system coordinates the depolarization of the ventricles, the 
QRS complex tends to look 'spiked' rather than rounded due to the increase in conduction velocity. A 
normal QRS complex is 0.06 to 0.10 sec (60 to 100 ms) in duration. Not every QRS complex contains a 
Q wave, an R wave, and an S wave. By convention, any combination of these waves can be referred to as 
a QRS complex. The duration, amplitude, and morphology of the QRS complex is useful in diagnosing 
cardiac arrhythmias, conduction abnormalities, ventricular hypertrophy, myocardial infarction, electrolyte 
derangements, and other disease states. Q waves can be normal (physiological) or pathological. The ST 
segment connects the QRS complex and the T wave and has a duration of 0.08 to 0.12 s (80 to 120 ms). It 
starts at the J point (junction between the QRS complex and ST segment) and ends at the beginning of 
the T wave. However, since it is usually difficult to determine exactly where the ST segment ends and the 
T wave begins, the relationship between the ST segment and T wave should be examined together. The 



Bioelectromagnetic activity of skeletal muscles is measured by non-invasive electromyo- 
graphy (EMG), usually recorded from surface electrodes placed on the skin above working 
muscleso An electromyograph detects the electrical potential generated by muscle cells, 
both when they contract or are at rest|^ A surface electrode may be used to monitor the 
general picture of muscle activation (see, e.g. [H]), as opposed to the activity of only a few 
fibres as observed using a needle. This technique is used in a number of settings; for ex- 
ample, in the physiotherapy clinic, muscle activation is monitored using surface EMG and 
patients have an auditory or visual stimulus to help them know when they are activating 
the muscle (biofeedback) o EMG signals are essentially made up of superimposed motor 
unit action potentials (MUAPs) from several motor units. For a thorough analysis, the 
measured EMG signals can be decomposed into their constituent MUAPs. MUAPs from 
different motor units tend to have different characteristic shapes, while MUAPs recorded 
by the same electrode from the same motor unit are typically similar!^ 

Bioelectromagnetism should not to be confused with the similar word bioelectromag- 
netics, which deals with the effect on life from external electromagnetism (see, e.g. |16]). 
Supported by the Bioelectromagnetics Society (with Dr. Stefan Engstrom as President) 
and the scientific journal with the same name, bioelectromagnetics is the study of how 
electromagnetic fields interact with and influence biological processes. Common areas of 
investigation include the mechanism of animal migration and navigation using the ge- 



typical ST segment duration is usually around 0.08 s (80 ms). It should be essentially level with the PR 
and TP segment. The normal ST segment has a slight upward concavity. Flat, down-sloping, or depressed 
ST segments may indicate coronary ischemia. ST segment elevation may indicate myocardial infarction. 

In contrast, to perform invasive intramuscular EMG, a needle electrode is inserted through the skin 
into the muscle tissue. The insertional activity provides valuable information about the state of the muscle 
and its innervating nerve. 

The electrical source is the muscle membrane potential of about -70mV. Measured EMG potentials 
range between less than 50 V and up to 20 to 30 mV, depending on the muscle under observation. Typical 
repetition rate of muscle unit firing is about 720 Hz, depending on the size of the muscle (eye muscles 
versus seat (gluteal) muscles), previous axonal damage and other factors. Damage to motor units can be 
expected at ranges between 450 and 780 mV. 

^"^A motor unit is defined as one motor neuron and all of the muscle fibers it innervates. When a motor 
unit fires, the impulse (called an action potential) is carried down the motor neuron to the muscle. The 
area where the nerve contacts the muscle is called the neuromuscular junction, or the motor end plate. 
After the action potential is transmitted across the neuromuscular junction, an action potential is elicited 
in all of the innervated muscle fibers of that particular motor unit. The sum of all this electrical activity 
is known as a motor unit action potential (MUAP). This electrophysiologic activity from multiple motor 
units is the signal typically evaluated during an EMG. The composition of the motor unit, the number of 
muscle fibres per motor unit, the metabolic type of muscle fibres and many other factors affect the shape 
of the motor unit potentials in the myogram. Nerve conduction testing is also often done at the same time 
as an EMG in order to diagnose neurological diseases. 

^^Notably MUAP size and shape depend on where the electrode is located with respect to the fibers and 
so can appear to be different if the electrode moves position. EMG decomposition is non-trivial, although 
many methods have been proposed (see, e.g. [15]). 



omagnetic field, studying the potential effects of man-made sources of electromagnetic 
fields, such as those produced by the power distribution system and mobile phones, and 
developing novel therapies to treat various conditions. 

Most importantly, in recent years a growing number of people have begun to use mo- 
bile phone technology. This phenomenon has raised questions and doubts about possible 
effects on users' brains. The recent literature review presented in [T7] has focused on 
the human electrophysiological and neuro-metabolic effects of mobile phone (MP)- related 
electromagnetic fields (EMFs) published in the last 10 years. All relevant papers have 
been reported and, subsequently, a literature selection has been carried out by taking sev- 
eral criteria into account, such as: blind techniques, randomization or counter-balancing 
of conditions and subjects, detail of exposure characteristics and the statistical analyses 
used. As a result, only the studies meeting the selection criteria have been described, 
evaluated and discussed further. This review has provided a clear scenario of the most 
reliable experiments carried out over the last decade and to offer a critical point of view 
in their evaluation. Its conclusion is that MP-EMFs may influence normal physiology 
through changes in cortical excitability and that in future research particular care should 
be dedicated to both methodological and statistical control, the most relevant criteria in 
this research field. 

More generally, epidemiological studies of radio frequency (RF) exposures and human 
cancers include studies of military and civilian occupational groups, people who live near 
television and radio transmitters, and users of mobile phones. Many types of cancer have 
been assessed in [18], with particular attention given to leukemia and brain tumors. The 
presented epidemiological results fall short of the strength and consistency of evidence 
that is required to come to a conclusion that RF emissions are a cause of human cancer. 
Although the epidemiological evidence in total suggests no increased risk of cancer, the 
results cannot be unequivocally interpreted in terms of cause and effect. The results are 
inconsistent, and most studies are limited by lack of detail on actual exposures, short 
follow-up periods, and the limited ability to deal with other relevant factors. For these 
reasons, the studies are unable to confidently exclude any possibility of an increased risk 
of cancer. Further research to clarify the situation is justified. Priorities include further 
studies of leukemia in both adults and children, and of cranial tumors in relationship to 
mobile phone use. 

In a recent technical paper [19], a new theoretical treatment of ion resonance phenom- 
ena was formulated around the Lorentz force equation. This equation defines the motion 
of a particle of mass m and charge q moving with a mean time r in a joint electric field E 
and magnetic field B 

mv -I- mw JT = (/(E -|- v x B), (1) 

where overdot means time derivative and x is the vector cross-product. The authors derive 
an expression for the velocity of a damped ion with arbitrary q/m under the influence 



of the Lorentz force on the right-hand side of ([T]). Series solutions to the differential 
equations reveal transient responses as well as resonance-like terms. An important result 
in this field is that the expressions for ionic drift velocity include a somewhat similar 
Bessel function dependence as was previously obtained for the transition probability in 
parametric resonance. The authors have found an explicit effect due to damping, so that 
previous Bessel dependence now occurs as a subset of a more general solution, including 
not only the magnetic field AC/DC ratio as an independent variable, but also the ratio 
of the cyclotronic frequency Q to the applied AC frequency to. The authors hypothesize 
that the selectively enhanced drift velocity predicted in their model can explain ICR- 
like phenomena as resulting from increased interaction probabilities in the vicinity of ion 
channel gates. 

This paper proposes rigorous geometrical bioelectrodynaTnics, which we define to be the 
common field theory underpinning both bioelectromagnetism and bioelectromagnetics. 

2 Exterior Geometrical Machinery in 4D 

2.1 From Green to Stokes 

Recall that Green's Theorem in the region C in (a;,y)— plane M^ connects a line integral 
!fdc (°^^^ ^^^ boundary dC of C) with a double integral JJ^ over C (see e.g., [20j) 

Pdx + Qdy = [[ (^- ^) dxdy. 
ac JJc\dx dy J 

In other words, if we define two differential forms (integrands of ^q^ and JJr') as 

1-form : A = Pdx + Qdy, and 

fdQ dP\ 

2-form : dA = ( — -— dxdy, 

\ox ay J 

(where d denotes the exterior derivative that makes a (p + 1)— form out of a p— form, see 
next subsection), then we can rewrite Green's theorem as Stokes' theorem: 

A = / dA. 

ac Jc 

The integration domain C is in topology called a chain, and dC is a ID boundary of a 2D 
chain C. In general, the boundary of a boundary is zero (see |35tl36j). that is, d{dC) = 0, 
or formally 9^ = 0. 



2.2 Exterior derivative 

The exterior derivative d is a generalization of ordinary vector differential operators {grad, 
div and curl see e.g., [31]) that transforms p— forms u into {p + 1)— forms du (see next 
subsection), with the main property: dd = d^ = 0, so that in M'^ we have (see Figures [T] 
and [2]) 

• any scalar function f = f{x,y, z) is a 0-form; 

• the gradient d/ = a^ of any smooth function / is a l~form 

cj = d/ = --dx + —-dy + -^dz; 
ox ay oz 



the curl a = dcj of any smooth 1-form cj is a 2-form 

„ ^- dydz + [-^ ^- dzdx + ^ ^- dxdy; 

ay oz J \oz ox J \ ox Oy J 

if cj = d/ ^ a = dd/ = 0. 

• the divergence /3 = da of any smooth 2-form a is a 3-form 

f dA dB dC\ 

f3 = da = ( — h — h — — ) dxdydz; if o; = do; =^ /3 = ddci; = 0. 

\ ox oy oz J 

In general, for any two smooth functions / = /(x, y, z) and g = g{x, y, z), the exterior 
derivative d obeys the Leibniz rule 



d{fg)=gdf + fdg, 

and the chain rule: 

d{gif))=g'{f)df. 

2.3 Exterior p— forms in W^ 

In general, given a so-called 4D cofranie, that is a set of coordinate differentials {dx^} £ W^, 
we can define the space of all p— forms, denoted r2^(M'^), using the exterior derivative 
d : r2P(M^) — > 0^"'""^(M^) and Einstein's summation convention over repeated indices (e.g., 
Ai dx^ = J2i=o ^i dx''), we have: 

1— form - a generalization of the Green's 1-form Pdx + Qdy, 

A = Aidx^ Erji(M^). 

For example, in 4D electrodynamics, A represents electromagnetic (co)vector poten- 
tial. 
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Figure 1: Basis vectors and one- forms in Euclidean M^— space: (a) Translational case; and 
(b) Rotational case [21]. For the same geometry in R^, see [35]. 
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Figure 2: Fundamental two-form and its flux in M^: (a) Translational case; (b) Rotational 
case. In both cases the flux through the plane u A f is defined as J J^^^cdpidq^ and 
measured by the number of tubes crossed by the circulation oriented by n A w [2T]. For 
the same geometry in M^, see [35j. 
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2— form - generalizing the Green's 2-form {dxQ — dyP) dxdy (with dj = d/dx^), 

B = dA G i7^(M'^), with components 

B = -Bijdx^ Adx^ , or 
2 

B = djAi dx^ A dx*, so that 

JDij = ZOjAi = OiAj OjAi = ^ji- 

where A is the anticommutative exterior (or, 'wedge') product of two differential 
forms; given a p— form a. G $7*'(IR^) and a g— form (3 G il'^(]R^), their exterior product 
is a (p + g)— form a A P ^ i7^^'^(M^); e.g., if we have two 1-forms a = aidx^ , and 
b = 6jdx*, their wedge product a A b is a 2-form a given by 

Q; = a A b = aibj dx^ A dx^ = —aibj dx^ A dx^ = —h A a. 

The exterior product A is related to the exterior derivative d, by 

d(Q; A /3) = da A /3 + {-If a A d/3. 



3 form 



C = dB ( = ddA = 0) G ri^(M^), with components 

C = —.Cijk dx' A dx^ A dx'^, or 

C = dkB[ij] dx^ A dx* A dx^ , so that 

Cijk = — 69fc-B[jj], where Bjjjj is the skew-symmetric part of Bij. 

For example, in the 4D electrodynamics, B represents the field 2-form Faraday, or 
the Lienard-Wiechert 2-form (in the next section we will use the standard symbol 
F instead of B) satisfying the sourceless magnetic Maxwell's equation, 

Bianchi identity : dB = 0, in components dkBuj] = 0. 



4 form 



D = 


= dC ( = ddB = 0) G rj^(M^), 


with components 


D = 


= diC[ijk] dx^ A dx"- A dx^ A dx'', 


or 


D = 


= -77 Dijki dx' A dx^ A dx'' A dx'', 


so that 


ijkl 


= -249zq,,,]. 
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2.4 Stokes Theorem in M^ 

Generalization of the Green's theorem in the plane (and all other integral theorems from 
vector calculus) is the Stokes Theorem for the p— form ld, in the nD domain C (which is 
a p— chain with a (p — 1)— boundary dC, see next subsection) 

u = I du:. 

dC Jc 

In the 4D Euclidean space M^ we have the following three particular cases of the Stokes 
theorem, related to the subspaces C of M^: 
The 2D Stokes theorem: 



The 3D Stokes theorem: 



The 4D Stokes theorem: 



A= / B. 



B= / C. 



C= / D. 

ac* Jc* 

2.5 Exact and closed p— forms and p— chains 

Notation: we drop boldface letters from now on. In general, a p— form (3 is called closed 
if its exterior derivative is equal to zero, 

d/3 = 0. 

From this condition one can see that the closed form (the kernel of the exterior derivative 
operator d) is conserved quantity. Therefore, closed p— forms possess certain invariant 
properties, physically corresponding to the conservation laws (see e.g., |30]). 
Also, a p— form f3 that is an exterior derivative of some {p — 1)— form a, 

13 = do, 

is called exact (the image of the exterior derivative operator d). By Poincare Lemma, 
exact forms prove to be closed automatically, 

dfj = d{da) = 0. 
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Since d? = 0, every exact form is closed. The converse is only partially true (Poincare 
Lemma): every closed form is locally exact. This means that given a closed p— form on 
a smooth manifold M (see Appendix), a G QP(M) on an open set U C M, any point 
m £ U has a neighborhood on which there exists a {pk — 1)— form /? G Q'P~^{U) such that 
dp = a\u. 

The Poincare lemma is a generalization and unification of two well-known facts in 
vector calculus: 
(i) If curl -F = 0, then locally F = grad /; and (ii) If div F = 0, then locally F = curl G. 

Poincare lemma for contractible manifolds: Any closed form on a smoothly contractible 
manifold is exact. 

A cycle is a chain C G Cp{M) such that dC = 0. A boundary is a chain C such that 
C = dB, for any other chain B £ Cp{M). Similarly, a cocycle (i.e., a closed form) is a 
cochain uj such that duj = 0. A coboundary (i.e., an exact form) is a cochain w such that 
uj = dO, for any other cochain 6. All exact forms are closed (lo = d6 ^ duj = 0) and all 
boundaries are cycles (C = dB =^ 5C = 0). Converse is true only for smooth contractible 
manifolds, by Poincare Lemma. 

2.6 Topological duality of p— forms and p— chains 

Duality of p— forms and p— chains on a smooth manifold M is based on the 



Period 

ic 



: fLo:={C,u), 
Jc 



where C is a cycle, a; is a cocycle, while (C, w) = u;{C) is their inner product {C,uj) : 
QP{M) X Cp{M) ^ M, a bilinear non-degenerate functional (pairing) on M |31j . From 
Poincare Lemma, a closed p— form to is exact iff {C,uj) = 0. 

The fundamental topological duality is based on the Stokes theorem, 



UJ = duo or {dC, u>) = (C, duo) , 

ac Jc 

where dC is the boundary of the p— chain C oriented coherently with C on M. While the 
boundary operator 9 is a global operator, the coboundary operator d is local, and thus 
more suitable for applications. The main property of the exterior differential, 

dod = d =0 =^ dod = d =0, (and converse), 

can be easily proved using the Stokes' theorem as 

= {d'^C,uo) = {dC,dLo) = {C,d^uo) = 0. 
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2.7 The de Rham cochain complex 

In the Euclidean 3D space M^ we have the following de Rham cochain complex 

^ n^{R^) -^ n\R^) ^ n^{R^) ^ n^{R^) -^ o. 

grad curl div 

Using the closure property for the exterior differential in M^, d o d = d? = 0, we get the 
standard identities from vector calculus 

curl • grad = and div • curl = 0. 

As a duality, in M^ we have the following chain complex 

^ Co{R^)^Ci{R^)^C2{R^)^C3{R^) ^ 0, 
(with the closure property d o d = d^ = 0) which implies the following three boundaries: 

Ci^Co = d{Ci), C2^Ci = d{C2), C3^C2 = d{C3), 

where Co G Co is a 0-boundary (or, a point), Ci G Ci is a 1-boundary (or, a line), 
C2 G C2 is a 2-boundary (or, a surface), and C3 G C3 is a 3~boundary (or, a hypersurface) . 
Similarly, the de Rham complex implies the following three coboundaries: 

C^^C^ = d{C^), C^ ^C^ = d{C^), C^^C^ = diC^), 

where C^ G 0,^ is 0-form (or, a function), C^ G fi^ is a 1-form, C^ G il^ is a 2-form, and 
C^ G r^^ is a 3-form. 

In general, on a smooth nD manifold M we have the following de Rham cochain 
complex (33] 

^ n^{M) ^ n\M) ^ n'^{M) -^ n^iM) -^ ^ o"(m) ^ o, 

satisfying the closure property on M, do d = d"^ = 0. 

2.8 (Co)homology on a smooth manifold M 

The subspace of all closed p— forms (cocycles) on a smooth manifold M is the kernel Ker(d) 
of the de Rham d— homomorphism (see Figure [3|), denoted by Z'f{M) C ilP(M), and the 
sub-subspace of all exact p— forms (coboundaries) on M is the image Im.{d) of the de Rham 
homomorphism denoted by BP{M) C ZP(M). The quotient space 

HP(M) ■= ^'^^'^ = Ker(rf:l^P(M)^QP+HM)) 
^ ^ ■ BPM lm{d : J^P-i(M) ^ Op(M)) ' 
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exact forms p -forms closed forms 

^/Ker(rf) 




(p — 1) -forms 



(p + 1) -forms 



Figure 3: A small portion of the de Rliam cochain complex, showing a homomorphism of 
cohomology groups. 

is called the pth. cohomology group of a manifold M. It is a topological invariant of a 
manifold. Two p— cocycles a, (3 £ QP{M) are homologous, or belong to the same cohom,ology 
class [a] e HP{M), if they differ by a (p - l)-coboundary a - (3 = dO £ Qp^^^M). 

Similarly, the subspace of all p— cycles on a smooth manifold M is the kernel Ker(9) 
of the 9— homomorphism, denoted by Zp[M) C Cp{M), and the sub-subspace of all 
p— boundaries on M is the image Im(5) of the 9— homomorphism, denoted by Bp(M) C 
Cp{M). Two p— cycles Ci,C2 S Cp are homologous, if they differ by a (p — 1)— boundary 
Ci — C2 = dB G Cp-i{M). Then Ci and C2 belong to the same hom.ology class [C] £ 
Hp{M), which is defined as 



Hp{M) 



Zp{M) _ Ker(5 : CpiM) -^ Cp_i(M)) 
Bp{M) ~ Im(9 : Cp+i{M) ^ Cp{M)) ' 



where Zp is the vector space of cycles and Bp C Zp is the vector space of boundaries on 
M. 



2.9 A cochain complex spanning the space-time manifold 



Consider a small portion of the de Rham cochain complex of Figure [3] spanning a space- 
time 4-manifold M, 

aP{M) -^ QP+\M) 



^P-\M)^'-'^^'^'-' ""^ 
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As we have seen above, cohomology classifies topological spaces by comparing two sub- 
spaces of n,'^: (i) the space of p— cocycles, Z^{M) := Kerdp, and (ii) the space of 
p— coboundaries, B^{M) := Imdp^i. Thus, for the cochain complex of any space-time 
4-manifold we have, 

(f = =^ BP{M) C ZP{M), 

that is, every p— coboundary is a p— cocycle. Whether the converse of this statement is 
true, according to Poincare Lemma, depends on the particular topology of a space-time 
4-manifold. If every p— cocycle is a p— coboundary, so that B^ and Z^ are equal, then the 
cochain complex is exact at 0P(M). In topologically interesting regions of a space-time 
manifld M, exactness may fail, and we measure the failure of exactness by taking the pth 
cohomology group 

HP{M) := ZP{M)/BP{M). 

2.10 Hodge star, codifferential and Laplacian 

The Hodge star operator :*: : r2^(M) — > r2"~P(M) on a smooth manifold M is defined as 



a A*/3 = /?A*Q = (q,/3)^, ioi a, 13 £nP{M). 

The -k operator depends on the Riemannian metric on M and also on the orientation (re- 
versing orientation will change the sign). The volume form fj, is defined in local coordinates 
on M as 

^ = vol = •(1) = Jdet{gij) dx^ A ... A dx", 

and the total volume on M is given by 



vol(M) = / *(1). 
Jm 



For example, in the 4D electrodynamics, the dual 2-form Maxwell -kF satisfies the electric 
Maxwell equation with the source [35], 

Dual Bianchi identity : d*F = *J, 

where *J is the 3-form dual to the charge-current 1-form J. 

For any to p— forms a,/3 G QP(M) with compact support, we define the (bilinear and 
positive-definite) Hodge L^— inner product as 

(a,/3) = / (a,/3)*(l) = / aA*/3. (2) 

Jm Jm 

We can extend the product (•,•) to L'^{QP{M)); it remains bilinear and positive-definite, 
because as usual in the definition of L^, functions that differ only on a set of measure zero 
are identified. 
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Now, the Hodge dual to the exterior derivative d : Q^{M) -^ QP^^{M) on a smooth 
manifold M is the codijjerential 6, a linear map 5 : r2P(M) -^ Qp^^{M), which is a 
generalization of the divergence, defined by |32j 

(^ = *d*(-l)"(P+^)+\ so that (i = *5*(-in. 

Applied to any p— form to G QP{M), the codifferential 6 gives 

If u; = / is a 0— form, or function, then 5f = 0. 

The codifferential 6 can be coupled with the exterior derivative d to construct the 
Hodge Laplacian A : i}P{M) -^ Qp{M), a harmonic generalization of the Laplace-Beltrami 
differential operator, given b}0 

A = 6d + d6. 

The Hodge codifferential S satisfies the following set of rules: 

• 66 = 6'^ = 0, the same as dd = d^ = 0; 

• 6A = A6, the same as dA = Ad; 

• 6* = {-1)P+^ ^ d; ir6 = {-l)Pird; 

• d6-k = i^6d; -kd5 = 6d-k; • A = A * . 

3 Classical Gauge Electrodynamics in 4D 

Recall that a gauge theory is a theory that admits a symmetry with a local parameter. For 
example, in every quantum theory the global phase of the wave ■0— function is arbitrary 
and does not represent something physical. Consequently, the theory is invariant under a 
global change of phases (adding a constant to the phase of all wave functions, everywhere); 
this is a global symmetry. In quantum electrodynamics, the theory is also invariant under 
a local change of phase, that is, one may shift the phase of all wave functions so that the 
shift may be different at every point in space-time. This is a local symmetry. However, in 
order for a well-defined derivative operator to exist, one must introduce a new field, the 
gauge field, which also transforms in order for the local change of variables (the phase in 
our example) not to affect the derivative. In quantum electrodynamics this gauge field is 
the electromagnetic potential 1-form A, in components given by 

A = Andx^, such that ^new = ^oid + df, 



''The difference d — 5 = dr, is called the Dirac operator. Its square d% equals the Hodge Laplacian A. 
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(where / is an arbitrary scalar function) - leaves the electromagnetic field 2-forni F = dA 
unchanged. This change df of local gauge of variable A is termed gauge transformation. In 
quantum field theory the excitations of fields represent particles. The particle associated 
with excitations of the gauge field is the gauge boson. All the fundamental interactions in 
nature are described by gauge theories. In particular, in quantum electrodynamics, whose 
gauge transformation is a local change of phase, the gauge group is the circle group U{1) 
(consisting of all complex numbers with absolute value 1), and the gauge boson is the 
photon (see e.g., [37|). 

In the 4D space-time electrodynamics, the 1-form electric current density J has the 
components J^ = [p-,]) = ip,jx,jy,jz) (where p is the charge density), the 2-form Faraday 
F is given in components of electric field E and magnetic field B by 



/ 



1^1/ 
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while its dual 2-form Maxwell -kF has the following components 
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3.1 Maxwell's equations 

The gauge field in classical electrodynamics, in local coordinates given as an electromag- 
netic potential 1-form 

A = A^dx^ = A^dx^ + df, (/ = arbitrary scalar field), 
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is globally a connection on a C/(l)— bundlqlj (see Appendix). The corresponding electro- 
magnetic field, locally the 2--forni 



F 
F 



dA, 
1 



in components given by 



Fr,^ dx^" A dx" , with F, 



flU 



duA^ 



d^A, 



is globally the curvature of the connection ^^^1 under the gauge-covariant derivative, 



where i 



Dfi = d^- ieAfj,, 
-1 and e is the charge coupling constant, 
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^^Recall that in the 19th Century, Maxwell unified Faraday's electric and magnetic fields. Maxwell's 
theory led to Einstein's special relativity where this unification becomes a spin-off of the unification of 
space end time in the form of the Faraday tensor [35j 

F = E Adt + B, 

where F is electromagnetic 2— form on space-time, E is electric 1— form on space, and B is magnetic 
2— form on space. Gauge theory considers F as secondary object to a connection-potential 1— form A. 
This makes half of the Maxwell equations into tautologies [3S], i.e., 

F ^ dA ^ dF = : Bianchi identity, 

but does not imply the second half of Maxwell's equations, 

5F = —AnJ : dual Bianchi identity. 

To understand the deeper meaning of the connection-potential 1— form A, we can integrate it along a 

7 
path 7 in space-time, x ^ y. Classically, the integral J A represents an action for a charged point 

particle to move along the path 7. Quantum-mechanically, exp lij A) represents a phase (within the 
unitary Lie group U{1)) by which the particle's wave-function changes as it moves along the path 7, so A 
is a C/(l)— connection. 

In other words. Maxwell's equations can be formulated using complex line bundles, or principal bundles 
with fibre U{1). The connection V on the line bundle has a curvature F = V^ which is a 2-form that 
automatically satisfies dF — and can be interpreted as a field-strength. If the line bundle is trivial with 
fiat reference connection d, we can write V = d -I- ^4 and F — dA with A the 1-form composed of the 
electric potential and the magnetic vector potential. 

^^The only thing that matters here is the difference a between two paths ^^ and 72 m the action J A 
|38) . which is a 2-morphism (see e.g., [231124] ) 




'if a gauge transformation is given by 



tp I 



■ e i> 
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Classical electrodynamics is governed by the Maxwell equations, which in exterior 
formulation read 

dF = 0, 6F = —AttJ, or in components, 

where comma denotes the partial derivative and the 1-form of electric current J = J^dx^ 
is conserved, by the electrical continuity equation, 

6 J = 0, or in components, J^,^ = 0. 

The first, sourceless Maxwell equation, dF = 0, gives vector magnetostatics and mag- 
netodynamics, 

Magnetic Gauss' law : divB = 0, 

Faraday's law : dtB + curlE = 0. 

The second Maxwell equation with source, 6F = J, gives vector electrostatics and elec- 
trodynamics. 

Electric Gauss' law : div E = 47r/3, 

Ampere's law : S^E — curlB = — 47rj. 

Maxwell's equations are generally applied to macroscopic averages of the fields, which 
vary wildly on a microscopic scale in the vicinity of individual atoms, where they undergo 
quantum effects as well (see below). 

3.2 Electrodynamic action 

The standard Lagrangian for the free electromagnetic field, F = dA, is given by |231 I24 1 [25] 

£(A) = i(FA*F), 



and for the gauge potential 

then the gauge-covariant derivative, 

transforms as 

and Dfj,il) transforms as 



e 



Df_,^ dfj, - \eAf_, ~ i{df_,K) 
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with the corresponding action 

s{a) = ^Ifa*f. 

Using the Hodge L^— inner product ([2]), we can rewrite this action as 



S{A) = l{F,F). (3) 



3.3 2D electrodynamics 



One of the reasons that gauge electrodynamics in 2D is especially simple is that the 
electromagnetic field F = dA, being a 2-form, can be written simply as a scalar field 
if = (p{x, t) times the volume form j44j 

F = (p vol. 

This means that the only Hodge duals we need to consider are the trivial ones, 

•k vol = 1 and * 1 = vol, 

so the 2D Lagrangian for vacuum electrodynamics is given by 



1 1 2 

L = —{if vol) A • ((^ vol) = —if vol A 1 = — vol. 

4 Quantum Electrodynamics in 4D 

Quantum electrodynamics (QED) is a relativistic quantum field theory of electrodynamics. 
QED was developed by four Nobel Laureates, P. Dirac, R. Feynman, J. Schwinger and 
S. Tomonaga, and F. Dyson, between 1920s and 1950s. It describes some aspects of how 
electrons, positrons and photons interact. QED mathematically describes all phenomena 
involving electrically charged particles interacting by means of exchange of photons. 

4.1 Dirac QED 

The Dirac equation for a particle with mass m (in natural units) reads (see, e.g., [25] ) 

{[-f^d^-m)ip = 0, (^ = 0,1,2,3) (4) 
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where i = \/— T ip{x) is a 4-coniponent spinoio wave-function, the so-called Dirac spinor, 
while j^ are 4x4 Dirac 'j— matrices, 
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They obey the anticommutation relations 

{7^7'^}=7^'^ + 7V = 25'^^ 

where ^^jy is the metric tensor. 

Dirac's 7— matrices are conventionally derived as 



7 



cr'^ 
-a'= 



(A: = 1,2, 3) 



where a are Pauli a— matrices^ (a set of 2 x 2 complex Hermitian and unitary matrices), 
defined as 



(Ti = (Ta; 



1 

1 



0-2 =0-^ 



-i 
i 



^■3 = (Tz 



1 

-1 



obeying both the commutation and anticommutation relations 

[o"j, Oj] = 2i Eijk CTk, {<7i, aj} = 25ij ■ I, 
where Eijk is the Levi-Civita symbol, 5ij is the Kronecker delta, and I is the identity 



matrix. 



The most convenient definitions for the 2-spinors, like the Dirac spinor, are: 



02 = 



and x^ 



X 



■^^In quantum mechanics, each Pauli matrix represents an observable describing the spin of a spin i 
particle in the three spatial directions. Also, iaj are the generators of rotation acting on non-relativistic 
particles with spin i . The state of the particles are represented as two-component spinors. 

In quantum information, single-qubit quantum gates are 2x2 unitary matrices. The Pauli matrices are 
some of the most important single-qubit operations. 
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Now, the Lorentz-invariant form of the Dirac equation (j4|) for an electron with a 
charge e and mass nie, moving with a 4-momentum 1-form p = p^dx^ in a classical 
electromagnetic field defined by 1-form A = A^^dx^, reads (see, e.g., }39t[25]): 

{ij^" [p^ - eAf,] - me} ^(x) = 0, (5) 

and is called the covariant Dirac equation. 

The formal QED Lagrangian (density) includes three terms, 

C{x) = Cem{x) + Cint{x) + Cc-pix), (6) 

related respectively to the free electromagnetic field 2-form F = F^j^dx^ A dx" , the 
electron-positron field (in the presence of the external vector potential 1-form A^^^), 
and the interaction field (dependent on the charge-current 1-form J = Jf^dx^). The free 
electromagnetic field Lagrangian in ([6]) has the standard electrodynamic form 

where the electromagnetic fields are expressible in terms of components of the potential 
1-form 

A = A^dx^^ by 

F^, = d^Al^' - d,A'°\ with 4°t = A^-t + A^. 

The electron-positron field Lagrangian is given by Dirac's equation ^ as 

£e-p(x) = ^{x) {ir [P^ - eAf] - me} ^{x), 

where V'(x) is the Dirac adjoint spinor wave function. 
The interaction field Lagrangian 

accounts for the interaction between the uncoupled electrons and the radiation field. 
The field equations deduced from ([6]) read 

{i7'^[p;.-e^-*]-me}V(^) = l^i'{x)A^, 

d^F^, = J,. (7) 

The formal QED requires the solution of the system ([7]) when A^(x), %l^{x) and ^(x) are 
quantized fields. 
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4.2 Feynman QED 

In Feynman's form of quantum electrodynamics |4U[|24[[^ . the process of quantization is 
transparent - it is performed using his path integral formalism, resulting in the quantum- 
field transition given (in natural units) by: 

{outAA\inA)= fv[A]e'^^^^ Wick f V[A]e~^^^\ (8) 

with action S{A) given by (J13p in the exponent (multiplied by the imaginary unit i). The 
Lebesgue integration in ([8]) is performed over electromagnetic fields A^, with the Lebesgue 
measure 



V[A]= \hn l[dA'^ (;u = l,..,4). 

The transformation Wick denotes the so-called Wick-rotation of the time variable t to 

> 

imaginary values t i-^ r = it, thereby transforming the complex- valued path integral into 
the real-valued (Euclidean) one. The absolute square of the transition amplitude ([8]) gives 
the transition probability density function, P{A) = \{outA\inA)\'^- Full discretization of 
([8|) ultimately gives the standard thermodynamic partition function 

Z = ^e-^V^, (9) 

j 

where E^ is the energy eigenvalue, T is the temperature environmental control parameter, 
and the sum runs over all j energy eigenstates. Prom @, we can further calculate all 
thermodynamical and statistical properties [H], as for example, transition entropy S = 
ks InZ, etc. 

The standard QED considers the path integral 

Z[A] = / T^IA] e' ^ ', where the action for the free e.-m. field is 

S[A] = Jd'x -\iF^uf =\j d^xA^{x){d''g^'''-d^'d^)A,{x). 

Z[A] is the path integral over each of the four spacetime components: 

V[A] = V[AfV[A]^V[AfV[Af. 

This functional integral is badly divergent, due to gauge invariance. Recall that F^j^^, and 
hence L, is invariant under a general gauge transformation of the form 

Af,{x) -^ A^{x) + -d^a{x). 
e 
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The troublesome modes are those for which Afj_{x) = d^a{x), that is, those that are gauge- 
equivalent to Af^{x) = 0. The path integral is badly defined because we are redundantly 
integrating over a continuous infinity of physically equivalent field configurations. To fix 
this problem, we would like to isolate the interesting part of the path integral, which counts 
each physical configuration only once. This can be accomplished using the Faddeev-Popov 
trick [42j, which effectively adds a term to the system Lagrangian and after which we get 



Z[A]= / V[A] exp 



d^x 



-T 



C-^^id'^A,) 



where ^ is any finite constant. 

Now, the observable O in the quantum field theory is a real-valued function of the 
gauge field: O : ^ ^ M, where A is the space of connections, so we can try to calculate its 
expected value using the electrodynamic action S{A) ([13]) as 



{n\TO{A)\n) 



J^V[A]0{A)e 

J^V[A]e-SiA) 



-S{A) 



In certain cases, it might be difficult to define exactly what this means. Even in the finite 
dimensional case, where 2? [A] is the Lebesgue measure on A, when we try to calculate (O) 
by the above formula we usually run into serious problems. To see this, let us do a sample 
calculation of the integral in the denominator, which is the partition function 



:= [ V[A]e- 
Ja 



S{A) 



(10) 



against which all other expected values are normalized. 

Thus, the Faddeev-Popov trick needs also to be applied to the formula for the two- 
point correlation function 



{n\TO{A)\n) 



lim 

T^oo(l-ie) 



fV[A]0{A)exp {j^r^d^xC 



/P[^]exp iJ'^^d^xC 
which after Faddeev-Popov procedure becomes 



{^\TO{A)\^) 



lim 

r-+oo(l-ie) 



JV[A]0{A)exY> i/^T^^a; C-^Adt'A^ 



J V[A] exp i /5^ d^x £ - ^ {^^'A 



The space of gauge transformations, ^, is a group which acts on the space A of con- 
nections by 

GxA^A, {(t>,A) h^ A + dcj). 
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A common way of eliminating divergences caused by gauge freedom is gauge fixing 
which means choosing some method to pick one representative of each gauge-equivalence 
class [A] and doing path integrals over these. There are problems with this approach, 
however. In a general gauge theory, it might not even be possible to fix a smooth, global 
gauge over which to integrate. But even when we can do this, the arbitrary choice involved 
in fixing a gauge is undesirable, philosophically. Gauge fixing amounts to pretending a 
quotient space of A, the space of connections modulo gauge transformations, is a subspace. 
A better approach is to use the quotient space directly. Namely, modding out by the 
action of ^ on ^ gives the quotient space A/Q consisting of gauge-equivalence classes of 
connections, and we do path integrals like 

Z=[ d[^]e-^([^]). 
JA/g 

When gauge fixing works, path integrals over A/G give the same results as gauge fixing. 
But integrating over A/G is more general and involves no arbitrary choices. 

Since gauge equivalent connections are regarded as physically equivalent, the quotient 
space A/Q of connections mod gauge transformations is sometimes called the physical 
configuration space for vacuum electrodynamics. A physical observable is then any real- 
valued function on the physical configuration space, O : A/Q -^ M, or equivalently, any 
gauge-invariant function on the space A of connections. In the case of noncompact gauge 
group even factoring out all of the gauge freedom may be insufficient to regularize our path 
integrals. In particular, there are certain topological conditions our space-time must meet 
for this programme to give convergent path integrals. As a consequence, in the case where 
the gauge group is M, we must sometimes take more drastic measures in order to extract 
meaningful results, and this leads to some interesting differences between the cases G = M 
and G = U{1) as the gauge group for electrodynamics. These differences are related to 
the famous 'Bohm-Aharonov effect' 



4.3 Cohomology criterion for path— integral convergence 



To eliminate the divergences in the path integral (jlOh it is necessary that we eliminate 
any gauge freedom from the space of connections A. In gauge theory, the action S{A) is 
invariant under the group of gauge transformations Q. In case of a gauge group ^ = R, 
we have gauge transformations of the form 

Ah^ A + dip, 

where 99 is an arbitrary (p — 1)— cochain. In other words, two p— connections are gauge 
equivalent if they differ by a p— coboundary dip G B^. Before we saw that the existence 
of gauge freedom implied path integrals such as the one above must diverge [43j. The 
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standard procedure of eliminating any gauge freedom from the space of connections A is 
to pass to the quotient space 

A/g = CP/BP, 

that is, p— connections modulo gauge transformations, effectively declaring gauge equiva- 
lent p— connections to be equal. In particuclar, the criterion in R electrodynamics for path 
integrals over the physical configuration space A/Q to converge is that the first cohomology 
be trivial [33]: 

\ converges for M electrodynamics J 

5 Electro— Muscular Stimulation 

In this section we develop covariant biophysics of electro-muscular stimulation, as an 
externally induced generator of the covariant muscular force, 

Fi = mija\ (11) 

where Fi are covariant components of the muscular force (torque) 1-form F = Fi{t)dx^, 
rriij are the covariant components of the musculo-skeletal metric/inertia tensor m = 
rriij dx^®dx^ , while a* = a*(t) d^i are contravariant components of the angular acceleration 
vector (see \12\ [23] for technical details). The so-called functional electrical stimulation 
(FES) of human skeletal muscles is used in rehabilitation and in medical orthotics to ex- 
ternally stimulate the muscles with damaged neural control. However, the repetitive use 
of electro-muscular stimulation, besides functional, causes also structural changes in the 
stimulated muscles, giving the physiological effect of muscular training. 

The use of low and very low frequency impulses in the body, delivered through elec- 
trodes, is known as transcutaneous stimulation of the nerves, electro-acupuncture and 
electro-stimulation. Here, an electromagnetic field accompanies the passage of the electric 
current through the conductive wire. This is generally known as the term 'electromagnetic 
therapy' H 



^^In the original sense acupuncture meant the inserting of needles in specific regions of the body. Electro- 
acupuncture supplies the body with low-volt impulses through the medium of surface electrodes to specific 
body regions or by non specific electrodes. Transcutaneous electric stimulation of the nerves (TENS) has 
for years been a well known procedure in conventional medicine. The impulses that are produced with this 
type of stimulation, are almost identical with those of electro-stimulation, yet many doctors still assume, 
that they are two different therapies. This has resulted in TENS being considered as a daily therapy, while 
electro-acupuncture or electro-stimulation were treated as 'alternative therapy'. Apart from the fact, that 
electro-acupuncture electro-impulses are delivered through needles, both therapies should be considered 
identical. Patients, who have reservations about the use of needles, can by the use of electric impulses 
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Figure 4: Schematic of electro-muscular stimulation {EMS). 

In particular, electro-muscular stimulation {EMS, see Figure S]) represents a union 
of external electrical stimulation fields, internal myofibrillar excitation-contraction paths. 



over surface electrodes on the skin, have a satisfactory alternative (see Figure |4]). We choose the term 
electro-stimulation, because the acupuncture system is not included in all therapies. Clinical tests have 
showed that there are two specific types of reactions: 

• The first reaction is spontaneous and dependent on the choice of body region. The stimulation of 
this part of the body results in an unloading, that can be compared with that of a battery. Normally 
this goes hand in hand with an immediate improvement in the patient. This effect of unloading 
may also be reached by non-specific electric stimulation. 

• The second normal reaction is of a delayed nature, that results in relaxation and control of pain. 
Moreover two other important effects follow, that begin between 10 and 20 minutes after the start 
of the treatment. This reaction is associated (combined) with different chemicals, such as beta- 
endorphins and 5-hydrocytryptamins. When using low and very low frequency stimulations, the 
second effect is obtained by the utilization of specific frequencies on the body. This is independent of 
the choice of a specific part of the body, because the connected electromagnet makes the induction 
of secondary electric current in the whole body possible. 
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and dissipative skin &: fat geometries, formally given by 

EMS = EMSf^Ms U EMSpaths \J EMSgeom- (12) 



In Feynman's QED formulation, corresponding to each of the three -BMS-phases in ()12p 
we formulate: 

1. The least action principle, to model a unique, external-anatomical, predictive and 
smooth, macroscopic EMS field-path-geometry; and 

2. Associated Feynman path integral, to model an ensemble of rapidly and stochastically 
fluctuating, internal, microscopic, fields-paths-geometries of the cellular EMS, to 
which the external-anatomical macro-level represents both time and ensemble av- 
eragec3 

In the proposed formalism, muscular excitation-contraction paths x*(t) are caused by 
electromagnetic stimulation fields $ = $ (t, x, y, z), by the the Lorenz force equation ([T]), 
while they are both affected by dissipative and noisy skin & fat shapes and curvatures, 
defined by the local Riemannian musculo-skeletal metric tensor mij (jlip . In particular, 
the electromagnetic 2-form Faraday F = Fa/sdxadx^, defines the Lorentz force 1-form 
Q = Qa dxa of the electro-muscular stimulation, 

Qa=PaOC qFa/3v'^, 

where q is total electric charge and v^ is the velocity vector-field of the stimulation flow. 
This equation says that the muscular force 1-form Qq generated by the simulation is 
proportional to the stimulation field strength Faf^, velocity of the stimulation flow v^ 
through the skin-fat-muscle tissue, as well as the total stimulation charge q. 

In the following text, we first formulate the global model for the EMS, to set up the 
general formalism to be specialized subsequently for each of the three EMS-phases. 

5.1 Global macro EMS-\e\e\ 

In general, at the macroscopic EMS-level we first formulate the total action S[^], our 
central quantity, which can be described through physical dimensions of Energy x Time = 
Effort, which is also the dimension of the Planck constant h {= 1 in normal units). This 
total action quantity has immediate biophysical ramifications: the greater the action - the 
higher the stimulation effect on the new shape. The action S'[$] depends on macroscopic 
fields, paths and geometries, commonly denoted by an abstract field symbol <&*. The 



^Recall that ergodic hypothesis equates time average with ensemble average. 
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EMS : A 



Figure 5: Simplified Feynman-like experimental approach to electrical muscular stimula- 
tion: the flow of electric current from the positive surface pad A® to the negative pad -B® 
can be approximated by the vector sum of complex vectors p^ e * , where 9^ are propor- 
tional to the time taken by each vector p^; this vector sum will be further developed into 
Feynman integral. 



action S\^] is formally defined as a temporal integral from the initial time instant tj^j to 
the final time instant t/jn, 

Sm = / £[$] dt, (13) 

with Lagrangian density, given by 

2.m = j d^xc{^\d,,^'), 

where the integral is taken over all n coordinates x^ = x^ (t) of the EMS, and d^j *^* are 
time and space partial derivatives of the <!>*— variables over coordinates. 

Second, we formulate the least action principle as a minimal variation 5 of the action 

dS[<^] = 0, (14) 

which, using variational Euler-Lagrangian equations, derives field-motion-geometry of the 
unique and smooth E'MS— transition map (or, more appropriately, 'functor' ^24j) 

T : STIMULf. => CONTRACTt ., => SHAPEt, , 
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acting at a macro-level from some initial time ij„j to the final time t/j„ (via the interme- 
diate time tmid)o 

In this way, we get macro-objects in the global EMS: a single electrodynamic stimu- 
lation field described by Maxwell field equations, a single muscular excitation-contraction 
path described by Lagrangian equation of motion, and a single Riemannian skin & fat 
geometry. 

5.2 Local Micro EMS-\e\e\ 

After having properly defined macro-level EMS, we move down to the microscopic cellular 
EMS-level of rapidly fluctuating electrodynamic fields, sarcomere-contraction paths and 
coarse-grained, fractal muscle-fat geometry, where we cannot define a unique and smooth 
field-path-geometry. The most we can do at this level of fluctuating noisy uncertainty, 
is to formulate an adaptive path integral and calculate overall probability amplitudes for 
ensembles of local transitions from negative EMS-pad A"^ to the positive pad B® (see 
Figure [5]). This probabilistic transition micro-dynamics is given by a multi field-path- 
geometry, defining the microscopic transition amplitude corresponding to the macroscopic 
-EMS'— transition map T. So, what is externally the transition map, internally is the 
transition amplitude. The absolute square of the transition amplitude is the transition 
probability. 

Now, the total E^MS"— transition amplitude, from the initial state STIMUL, to the 
final state SHAPE, is defined on EMS 

{SHAPE\STIMUL) : STIMULt^ ^ SHAPEt,, (15) 

given by modern adaptive generalization of the Feynman's path integral, see [25]). The 
transition map (J15p calculates overall probability amplitude along a multitude of wildly fluc- 
tuating fields, paths and geometries, performing the microscopic transition from the micro- 
state STIMULtQ occurring at initial micro-time instant to to the micro-state SHAPEt^ 
at some later micro-time instant ti, such that all micro-time instants fit inside the global 



■^^Here, we have in place n— categorical Lagrangian-field structure on the muscular Riemannian config- 
uration manifold M [24] 

using 



with 



$* 


:[0,1]^M, $':$^^$1, 




, ^ ^ f dc \ dc 



[a;o,a;i] ^ [<&o,*'i] 
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transition interval tQ,ti, ...,ts G [tini,tfin]- It is symbolically written as 

{SHAPE\STIMUL)= I V[w^e'^^'^\ (16) 

where the Lebesgue integration is performed over all continuous $co„ = fields + paths + 
geometries^ while summation is performed over all discrete processes and regional topolo- 
gies $;^j^. The symbolic differential 2?[u;<5] in the general path integral (J16p . represents an 
adaptive path measure, defined as a weighted product 



'D[w^\= lim TTty^d^*, {i = l,...,n = con + dis), (17) 

N — >oo -'--'- 

which is in practice satisfied with a large A^. 

In this way, we get a range of micro-objects in the local EMS at the short time-level: 
ensembles of rapidly fiuctuating, noisy and crossing electrical stimulation fields, myofib- 
rillar contraction paths and local skin & fat shape-geometries. However, by averaging 
process, both in time and along ensembles of fields, paths and geometries, we can recover 
the corresponding global, smooth and fully predictive, external EMS transition-dynamics 

r. 

5.3 Micro— Level Adaptation and Muscular Training 



The adaptive path integral (|16H17p incorporates the local muscular training process ac- 
cording to the basic learning formula 

NEW VALUE = OLD VALUE + INNOVATION, 

where the term VALUE represents respectively biological images of the STIMUL, 
CONTRACT and SHAPE. 

The general synaptic weights Wg = Wsit) in pT|) are updated by the homeostatic neuro- 
muscular feedbacks during the transition process T, according to one of the two standard 
neural training schemes, in which the micro-time level is traversed in discrete steps, i.e., 
if i = to,ti, ...,ts then f + 1 = ti,t2, •••,*s+i: 

1. A self-organized, unsupervised, e.g., Hebbian-like training rule (see, e.g. [28j: 

Wsit + 1) = Wsit) + -{w^sit) - <(i)), (18) 

rj 

where a = cr{t), 77 = r]{t) denote signal and noise, respectively, while superscripts d 
and a denote desired and achieved muscular micro-states, respectively; or 
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2. A certain form of a supervised gradient descent training: 

Wsit + l) = Ws{t)-vVJ{t), 



(19) 



where 77 is a small constant, called the step size, or the training rate, and V J{n) 
denotes the gradient of the 'performance hyper-surface' at the i— th iteration. 

For more details EMS-fields, paths and geometries, see 



6 Appendix: Manifolds and Bundles 

6.1 Manifolds 

Smooth manifold is a curved nD space which is locally equivalent to M". To sketch it 
formal definition, consider a set M (see Figure [6]) which is a candidate for a manifold. 
Any point x G iWo has its Euclidean chart, given by a 1-1 and onto map ip^ : M ^ W^, 
with its Euclidean image Vi = (pi{Ui). More precisely, a chart c/9j is defined by 

iPi-.MDUiBx^ iPiix) eViC M."-, 

where Ui C M and Vi C M" are open sets. 




Figure 6: Geometric picture of the manifold concept. 

Clearly, any point x £ M can have several different charts (see Figure [6|). Consider a 
case of two charts, if^, ipj : M ^ R", having in their images two open sets, Vij = ipi{Uir\Uj) 
and Vji = (Pj{Ui D Uj). Then we have transition functions (p^j between them, 

^1 



^ij = V>j ° ^i 



Vi 



V 



V 



Jll 



locally given by p>ij{x) = (fjiPi ^{x)). 



^^Note that sometimes we will denote the point in a manifold M by m, and sometimes by x (thus 
implicitly assuming the existence of coordinates x = (k')). 
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If transition functions ip^j exist, then we say that two charts, (p^ and 99 • are compatible. 
Transition functions represent a general (nonhnear) transformations of coordinates, which 
are the core of classical tensor calculus. 

A set of compatible charts ^p^:M^ M", such that each point x £ M has its Euclidean 
image in at least one chart, is called an atlas. Two atlases are equivalent iff all their charts 
are compatible (i.e., transition functions exist between them), so their union is also an 
atlas. A manifold structure is a class of equivalent atlases. 

Finally, as charts ip^ : M ^ M" were supposed to be 1-1 and onto maps, they can be 
either homeomorphisms, in which case we have a topological (C^) manifold, or dijjeomor- 
phisms, in which case we have a smooth {C^) manifold. 



6.2 Fibre Bundles 

On the other hand, the well-known tangent bundle TM and cotangent bundle T*M 
of a smooth manifold M, respectively endowed with Riemannian geometry (suitable for 
Lagrangian dynamics) and symplectic geometry (suitable for Hamiltonian dynamics) - 
are special cases of a more general geometrical object called fibre bundle. Here the word 
fiber F of a map tt : Y ^ X denotes the preimage 7r~^(a;) of an element x £ X. It is 
a space which locally looks like a product of two spaces (similarly as a manifold locally 
looks like Euclidean space), but may possess a different global structure. To get a visual 
intuition behind this fundamental geometrical concept, we can say that a fibre bundle Y 
is a homeomorphic generalization of a product space X x V (see Figure [7]), where X and 

V are called the base and the fibre, respectively, tt : F — > X is called the projection, 
Yx = n~^ (x) denotes a fibre over a point x of the base X, while the map / = vr^^ : X ^ Y 
defines the cross-section, producing the graph {x, f{x)) in the bundle Y (e.g., in case of a 
tangent bundle, f = x represents a velocity vector-field, so that the graph in a the bundle 

Y reads {x,x)). 



bundle Y^XxV 




V 



Y 


l^ 


f] 


\7t 



X 



X 



Figure 7: A sketch of a fibre bundle Y Ki X x V as a generalization of a product space 
X X V; left - main components; right - a few details (see text for explanation). 
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The main reason why we need to study fibre bundles is that all dynamical objects 
(including vectors, tensors, differential forms and gauge potentials) are their cross-sections, 
representing generalizations of graphs of continuous functions. 
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